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Abstract. We consider solutions of the stochastic equation X = X)j=i AiXi + B, where N is a 
random natural number, B and Ai are random positive numbers and X; arc independent copies 
of X, which arc independent also of N,B,Ai. Properties of solutions of this equation arc mainly 
coded in the function m(s) = E [ J^—j Af]. In this paper we study the critical case when the 
function m is tangent to the line y = 1. Then, under a number of further assumptions, we prove 
existence of solutions and describe their asymptotic behavior. 



1. Introduction 

The main purpose of the present paper is to study a class of linear stochastic equations, existence 
of their solutions and to describe properties of those solutions. The simplest example we have in mind 
is the random difference equation, called often also a hrst order random coefficients autorcgressive 
model, 

(1.1) X = AX + B, 

where all the random variables are real valued, X is independent of the pair (A, B), and the sign 

'=' denotes equality in distribution. It is well-known that the equation above has a unique solution 
if E[ log A] < and E[log + < oo. The celebrated Kestcn theorem [19] says that if a random 
variable X is a solution of the equation (1.1), then under a number of assumptions, the main being 
existence of a positive a such that E[A Q ] = 1, the random variable X is a-regularly varying, i.e. 

lim x a P[X > x] = C+, 

x — > oo 

for some positive constant C + (see also the paper of Goldie [15] for a transparent and elegant proof). 
Since the random difference equation appears both in many applied models e.g. in financial mathe- 
matics and in purely mathematical problems, the last result found enormous number of applications 
in the literature. 

In this paper we consider general linear stochastic equations, i.e. equations of the form 

N 

(1.2) X = ^AiXi + B, 

i=l 

where X, X^ are i.i.d. and independent of (N, B,Ai,A 2 ---). Notice that the last formula depends 
only on TV first values of Aj's, therefore without any loss of generality, we assume that Ai = for 
i > N. Moreover in this paper we restrict our attention to positive random variables, i.e. we assume 
that Xi, Ai and B are positive. 
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Equation (1.2) is also called the inhomogeneous smoothing transform and the explanation of 
this name is the following. Given [i 6 'P(R + ) we define Tfi as the law of X^=i AiXi + B, where 
Xi,X 2 , ... is i.i.d. sequence with distribution fi, independent of the vector (N, B, A\, A 2 , . . . ). Then 
any fixed point of T is characterized as a distribution of a random variable X that satisfies (1.2). 

This equation has gained importance in the last few years, since it turns out to be closely related 
to important objects in the computer science: the Quicksort algorithm [21, 24] (and other divide 
and conquer algorithms), the Pagerank algorithm [26, 17, 18] (being in the heart of the Google 
engine) and in stochastic geometry [22] . The inhomogeneous equation was recently used to describe 
equilibrium distribution of a class of kinetic models see e.g. [7] 

Equation (1.2) is also a generalization of the homogeneous smoothing transform, which is defined 
exactly as above but without the inhomogeneous term B 7 i.e. with B = a.s. Thus, we say that X 
is a solution (or a fixed point) of a homogeneous smoothing transform if 

N 

(1.3) I = 5]Mi, 

i=i 

where Xi, X 2 , ■ ■ ■ are independent copies of X and the vector (N, Ai,A 2 ,..) is independent of the 
sequence {-?Q}. The last equation appeared in the literature already in the eighties in connection 
with studying interacting particle system [13]. It turned out also to have a number of applications 
e.g. in branching random walks [16]. 

Existence of solutions of (1.3) and their properties were deeply studied in [13, 20] (see also the 
recent paper [3]) and it turns out that their properties are encoded in the function 



(1.4) m(t) = E 



N 

£4 



Notice, equation (1.3) does not have a unique solution since tX for t G R also solves it as long as X 
does. We summarize known results in the following Lemma 

Lemma 1.5 ([13, 20]). If EN > 1 and inf s e m i] m(s) < 1, then the set of solutions of (1.3) is 



< oo, E [./V 1+<5 ] < oo and for some a € (0, 1) we have 



nonempty. Moreover if E 
m(a) = 1, m'(a) < then 

lim P [X > x] x a = c if m'(a) < 0, 

x— *-oc 

lim P [X > x] x a l log x = c if to' (a) = 0, 

x— loo 

for some positive constant c. 

We begin the study of the nonhomogencous smoothing transform explaining how to construct a 
solution to equation (1.2) (see [5, 4] for more details). Let T = Ufe>o ^ ^ e an mnn itc Ulam-Harris 
tree, where N° = {0}. For v = (i 1; . . . , i n ) we define the length \v\ = n and by vi we denote the 
vertex 12, ■ ■ ■ ,i n ,i). We write u < v if u is a proper prefix of v, i.e. u = (i\, .., ik) for some k < n. 
Moreover we write u < v if u < v or u = v. Now we take {(B(v), A\(v), A 2 {v), . . . )} v gT a family 
of i.i.d. copies of (B, Ai, A 2i ■ ■ . ) indexed by the vertices of T. For v € T we also define a random 
variable L{0) = 1 and L(vi) = L(v)Ai(v). We can define now 

(1.6) R=J2 L (v)B(v). 

One can easily check that if the series above is finite almost surely then the random variable R 
satisfies (1.2). However also the converse is true. Alsmeyer and Meiners [5] proved that existence 
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of solutions of (1.2) is equivalent to finiteness of the series (1.6). Knowing that there exists one 
solution, one can construct a whole family of solutions just adding to R any Y being a solution 
of (1.3). However R is distinguished by the property that it is the minimal solution (in the sense 
of stochastic domination i.e. P [R > t] < P [X > t] for any other solution X), sec [5] for more 
details. Another useful property is that R is the unique solution that is measurable with respect 
to the input data (B(v),A\(v),A2(v), . . . ) ve j- (compare with [2] where it is called an endogeneous 
solution) . From now we call R the minimal solution. Therefore, if we can describe the tail of R, in 
view Lemma 1.5, wc obtain a full description of tails of all solutions (1.2). 

Similarly like in the homogeneous case the fundamental role in description of solutions of (1.2) 
plays the function m defined in (1.4). The necessary condition ensuring finiteness of (1.6) is that 
m{to) < 1 for some to € [0, 1]. However sufficient conditions are still not established. It is known 
[5, 17] that if m(s) < 1 for some s G (0, 1) and EB S < oo then R is well defined. 

Jelenkovic and Olvera-Cravioto [17, 18] proved that R has a power law distribution: 

Lemma 1.7. Let (B, Ai, A2, . . .) be a nonnegative random vector, with iVGNU{oo},P[.B>0] >0 
and R be the minimal solution to (1.2) given by (1.6). Suppose that 
• the equation m(s) = 1 has 2 solutions a < ft; 



< 00. 



E [B f3 ] <oo,0< m'(P) = E [Ef = i K log A, 
there exists j > 1 with P [N > j, Aj > 0] > such that th 



e measure i 



[log A, £ du, A 3 > 0, N > j] 



is nonarithmetic; 
In addition, assume that 



a) m(l) = 



Ei=i Ai 



< 1 and E 



(E*Li A 



< 00, ifP> 1; 



b) E 



0/(1- 



l+e 



< for some < e < 1, if < ft < 1. 



The 



[R > t] ~ Ct 



t —> 00, 



for some C > 0. 

Remark 1.8. Positivity of the constant C was not discussed in [17, 18] and was proved recently in 
[12]. 

Summarizing if a < ft are two solutions of the equation m(s) = 1 and a < 1, then the minimal 
solution R of (1.2) has a power law of index j3. Any other solution X of (1.2) satisfies P [X > t] ~ 
Ct~ a . 

The main purpose of the present paper is to complete the picture and to study the critical case, 
when the equation m(s) = 1 has exactly one solution a < 1 and then m'(a) = 0, i.e. when the 
graph of the function m(s) is tangent to the line y = 1. For the random difference equation (1.1) 
this corresponds to the situation when the graph of the Mellin transform EL4 S ] is tangent to the 
line y = 1 at 0, i.e. when E[logA] = 0. Then it is known that equation (1.1) has no solutions, 
nevertheless when written in terms of measures has solutions in the class of Radon measures on M. 
Existence and asymptotic properties of solutions were studied in [6, 8, 9]. For the homogeneous 
smoothing transform the critical case was considered by Durrett, Liggett [13] and Liu [20] and is a 
part of Lemma 1.5 (see also [10] for the case a = 1). 

The only result we know concerning the inhomogeneous smoothing transform in the critical case is 
due to Alsmeyer and Mciners [5], who proved that for a < 1/5 (and under some further assumptions) 
the series (1.6) is finite, providing thus a solution to (1.2). 

The main result of this paper is the following 
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(1.10) There exists < a < 1 such that m(a) = E A" = 1, 




(1.11) m'(o)=E ^Aflog^ =0, 



_i=i 



(1.12) E[iV] > 1, 

(1.13) For some j the measure P [log Aj £ du, Aj > 0,N > j] is nonarithmetic, 

T N 1 



(1.14) E 7V 1+lS + B a+S + (Ar 5 + A? +<5 ) < oo, for some < 5 < 1 -a. 



j=i 

TTierc t/ie minimal solution R of (1.2) is weH defined and moreover 
(1.15) lim £ Q P [i? > t] = C+ 



anc! t/ie constant C+ is strictly positive. 

Thus in the critical case the tail of the minimal solution of (1.2) is of the order t~ a , whereas the 
tails of all the other solutions behaves at infinity like logi t~ a . 

We finish the introduction with an overview over the organization of the paper. In Section 2 we 
prove that F[R > t] < Ct~ a , that in particular implies that R is finite a.s. The most essential part 
of the proof is contained in Section 3. We reduce the problem to study behavior at of the Laplace 
transform (j> of R. Considering as a solution of the Poisson equation we first prove that it behaves 
regularly at (Section 3.1) and the deduce the correct asymptotic (Section 3.3). Finally, applying 
some arguments based on the Landau theorem and holomorphic functions, we prove positivity of 
the limiting constant (Section 3.4). 

The authors are grateful to Jacek Zicnkiewicz for stimulating discussions on the subject of the 
paper. 



where R is the random variable defined in (1.6). In particular R is finite a.s. 

Corollary 2.2. E[R ] is finite for every (3 < a. 

We start with the following lemma 

Lemma 2.3. Let {Y n } n£ jq be a sequence of i.i.d. random variables such that EYi = 0. Let S n = 
Yli—i Yi be the sequence of the partial sums. Then, for any strictly positive constant S, the function 



2. Existence of a solution and its first estimates 



In this section we prove 



Proposition 2.1. Assume hypotheses of Theorem 1.9 are satisfied, then 

F[R >t}< Ct~ a , 




is bounded. 
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Proof. By definition the function W can be nonzero only for positive x. Let L = inf{i : Si < 0}, 
then 



W(x) = E 



E 



e L [S j +x>0 for j<i] 



E 



E 



i=0 
L-l 

6 ->-[S i +x>0 for j<i] 

i=0 
L-l 

E< 

i=0 
L-l 



■E 



6 J-[S 3 +a;>0 for j<i] 



-<5(x+Si) 



[x>0] 



E 



e L [Sj+x>0 for L<j<i 



i=0 



e-^l [x > 0] +E[VF(x + 5 L )] 



Notice that the first expression above is just a finite constant, since by the reflection principle [14] 

"L-l 



E 

i=0 



,-SSi 



E- 

n=0 



=: C < oo, 



where T„ is the sequence of upward ladder times: To = 0, T n = > T„_i : Si > S't„_ 1 }- 

Moreover the function f(x) = e~ Sx l[ x > ] is directly Riemann integrable (dRi), i.e. it is integrable 
and satisfies 

(2.4) ! im/i E su p \f( x )-f(v) l = °. 

where I n (/i) = (nh, (n + l)h]. 

Finally, we can express W as the convolution of the function / with the potential of the transient 
random walk V n , where V n is the sum of n independent copies of Sl- Therefore, independently of 
x, we have 



oo oo 



[x+v„e[k,k+i)] 



n=0 fc=0 



< e~ Sk ^f[x-k + V n e [0, 1)] <C'^2 e ~ Sk < oo. 

fc=0 11=0 fe=0 

where the uniform bound in the last line follows from Proposition 2.1 in [23, CH. 5] 
Let us introduce a random variable Y with distribution given by 



□ 



(2.5) 



E[/(F)] =E 



N 



J2f (-log AJA? 



for any positive Borel function /. By (1.10), the right hand side indeed defines a probability measure. 
The main properties of Y , we arc going to use, are summarized in the following lemma 

Lemma 2.6. The random variable Y is centered (EY = 0), nonarithmetic (the closed subgroup 
generated by the support of the measure P [Y G dx] is TBI) and has finite exponential moments 



E [e ±5Y ] < oo 



for some 6 > 0. 
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Proof. We have 



E [e ±SY ] = E 



N 



E A 

.1=1 



e [r] = e 



AT 



m(a =f 6) < oo, 



0. 



Nonaxithmecity follows from assumption (1.13). □ 

Proof of Proposition (2.1). We compare behavior of the sum R = X^eT L(v)B(v) with behavior of 
the maximum R — max„eT L(y). 
We first prove that 

(2.7) F[R > t] < cr a , 

for some positive constant C. 

Let {Yi} be a sequence of independent copies of Y defined in (2.5) and let S n be the sequence of 
their partial sums. Applying the definition of Y and reasoning by the induction (see e.g. [1]) one 
can easily prove the following well-known formula: 



(2.8) 



E[e Q5 "/(5i,...,5„)] = 



E /(- lo g-^( w i))->-logi(w„)) 

\v\—n 

valid for a fixed n and any test function / : K" — > M. 

Putting f(x!,..,x n ) = l[ a!l >_iogt,..., a! „_ 1 >-iogt,x n <-logt] we obtain 



R > t 



P [L(v) > t, for some v £ T) 

P > t for some |w| = n and £( u ) < t for w < y] 



< 



E E 



[— log L(vi)>— log t,...,— log L(ti„_i)> — log t,— log L{v)<— log t] 



— E^ [ e "' S " : '-[Si>-logt,...,S„-.i>-logt,S„<-logt]] 

n 

< t~ a , 

hence we obtain (2.7). 
Next we write 

P [R > t] < F[R > t] + P {R > t) n {R < t) . 
In view of (2.7) it is sufficient to estimate only the second term. Taking 7 = a + S < 1, we have 

v\{r >t}n{R<t} 



< 



< 



E L ( U )' B ( U ) 1 [i(f')<* for «'<*>] > * 
]T L( V )^( V ) 7 l[L (i/)<t for u'<t;] 



- tl L J 



E ^( w ) 7 l[L(«')<t for «'<«] 



«6T 
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Applying again (2.8) we obtain 



E 



^ ^(v) 7 l[L(t,')<t for v'<v) 



= 5Z E X] jy ( t, ) 7l [£(f')<* for v'<v 
n L |tj|=n 

= ^E[e aS ™e 7S "l[s fc +iogt>o for fc<n]] 
n 

n 

= ^PF(logt), 



L [S fc +logt>0 for fe<n] 



where IF is the function defined in Lemma 2.3, which as we already know is bounded. Finally, since 
E[£? 7 ] < oo, we obtain 



{R > t) n {r < t) 



< ct~ a . 



□ 



3. Tail of the solution 

3.L The Poisson equation. For a non-negative random variable X by <fix(t) = E [e _tx ] we 
denote its Laplace transform. For simplicity we write <f> = 4>r for the Laplace transform of R. To 
prove our main result we use the duality between the tail behaviour of R and the behaviour of its 
Laplace transform (f> near given by the following Tauberian theorem (its proof can be found e.g. 
in the book of Feller [14], Example c) after Theorem 4 in Chap. XIII). 

Lemma 3.1 (Tauberian Theorem). For < a < 1 and a slowly varying function L the following 
are equivalent: 



i) 
it) 



,. x a P\R>x] 
hm i— — = 1 



L(x) 

i - m 

tHMO t a L(\/t) ' 



lim 



r(i-a). 



Thus, in order to describe the tail of R, i.e. V[R > t], it is sufficient to study its Laplace transform 
4> and prove 

hml^ = C. 

t-yO t a 

It is convenient for our purpose to change the coordinates and define 

D(x) =e ax {l-(f){e- x )). 

Then our aim is to prove 

(3.2) lim D{x) = C. 

X— ¥00 

We will often use the following well-known lemma 

Lemma 3.3. For any positive random variable X and any < 7 < 1 we have 

i-<m*) <m- 7 )E[x^. 
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Proof. Notice that E [f(X)\ = J °° /'(s)P [X > s] ds for nonnegative X and any differentiable, mono- 
tone function / such that /(0) = 0. Then by Chebyshev inequality 

POO 

1 - <t>x (t) = E [l - er tx ] = (1 - e~ s )'P [tX > s] ds 

Jo 

/>oo 

< / e- s E[{tXy]s-" t ds^r(l- 1 )E[X' r }f f . 
Jo 

□ 

To prove (3.2) we apply some techniques described in the paper of Durrett and Liggett [13], who 
considered solutions of the homogeneous equation (1.3). We adopt their ideas, however it turns out 
that adding the additional term B causes serious problems, hence we will present here all the details 
of the proof. 

Let Y be the random variable defined in (2.5). We consider D as a solution of the Poisson 
equation. 

Lemma 3.4. The function D satisfies the following Poisson equation 

E[D(x + Y)] = D(x)+G(x), 

where 



G(x) = e ax E 



N 



N 



£(l-0( e -^))- 1-e 



e~ x B 



i=l 



i=l J 



Proof. Notice first that rewriting equation (1.2) in terms of Laplace transform cf> we obtain 



<p(t) = E 



N 



-ID 



Hence by the definition of D and the equation above we have 

E [D(x + Y)] - D{x) = E e a{x+Y) (l - <p (V (:r+Y) )) - e ax {l - <p{e- x )) 



e ax E 



' N / N \ " 

e-^ogAi ^_ (f) (g-a+iog*)) Af-(l- e~ e ~ XB H ^(e- x Ai) 

. i=l V i=l / . 

' N I N \ ' 

= e ax E 2(1 - <f>(e- x A t )) - 1 - e~^ B J] ^{eT'Ai) 

,i=l V i=l / . 



G(x) 



□ 



We need also the following technical lemma saying that for any t € (0, a + S) the sum X)i=i 
that appears under the expected value in the definition of m (1.4), has moment bigger than 1. 

Lemma 3.5. Assume that E [iV 1+(5 ] < oo for some S > 0. Let Xi be an arbitrary sequence of 
random variables. Then for any r > 1 and any pff (l, r ^+p ) we have 



E 



N \ Pi 



\i=l 



N 



where CV iP is a constant depending on r and p. 
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Proof. We will use the following well-known inequality (being just a simple consequence of the 
Holder inequality) 

(n \ r n 

i=i / i=i 

where a, > 0. 

Applying first the last inequality and then the Holder inequality with parameters r/p and (r/p)' 
(given q > 1, we denote by the conjugate real number such that 1/q + 1/q' = 1) wc obtain 







[{p'"f. 


< E 







E 



Notice that since 



A? 



p/r 



N {r-l)v/r I^Xi 



< E 



' (r-l)f fry 



l/(r/p)' . 



N 



E^ 

i=i 



/V' 



(r-1) 



p/r 



r \p 



<l + 5 



in view of our assumptions both expressions above are finite. □ 
We prove now a weaker result than (3.2) saying that the function D behaves regularly at infinity. 

Proposition 3.6. Assume that hypotheses of Theorem 1.9 are satisfied. Then for any y G R we 

have 

lim D t\ V) = 1. 

x->oo D[x) 

Proof. We divide the proof of the proposition into several steps. Assume first additionally that the 
following condition is satisfied 



(3.7) 



1 - E \e- tB ] 
lim l — — l - = 

4-V0 1 - <f>(t) 



Step 1. First wc will show that 
(3.8) 

We can write 



G(x) 



lim _ 

x^ao D(x) 



(3.9) 



lim E 

x— >oo 



lim E 



x^oo D{X) 



1) 



f - cj)(e- x ) 

nti ^(e" x A) - 1 + EL(1 - <^(e-^i)) 



f 



(e—) 



By our assumptions the first term is equal to since the Laplace transform (f> is bounded by 1. 
In order to show that the second limit (3.9) is zero we will use the following inequality valid for 

< Ui < Vi < 1 



N 



N 



N 



N 



(3.10) 



~[ Ui -l + ^(1 - Ui ) > Y[ v, - 1 + - Vi) (see [13, (2.5)] 



i=i 



Next we will deduce that the expression under the expectation in (3.9) is positive. 



10 



D. BURACZEWSKI, K. KOLESKO 



In order to bound this limit from above we use the inequality u < e ^ ") for uel. Therefore, 
we can write 



(3.11) 



E 



JV 



JV 



n0(e- x A 4 )-l + ^(l-0(e^)) 

N \ 

£(1 - <Ke-*Ai)) - 1 + £(1 - 0(e-^)) 



< E 



cxp 



i=i 



i=l 



Now we will split the problem into two separate cases 



Step 1, case i) There exists constant M such that for any i, Ai < M a.s. Since is a Laplace 
transform of a non-negative random variable, (1 — <p(u))/u is a decreasing function whereas (1 — 4>{ u )) 
is increasing. Hence 

1 - 4>(e-*Ai) < max(A i , 1)(1 - <KO), 

and thus 

AT JV 

2 [1 - ^(e-A-)] < (JV + J] 4)(1 " 0(e" x )). 

i=l i=l 

Therefore, since the function i^(u) = e~" — 1 + u is increasing on [0, oo), F(u)/u is bounded and 
tends to as u — >■ we can apply the Lebesgue theorem and obtain 



3 



lim sup ■ 



Utx tte-'At) 1 + Eili(l - ^(e- a ^)) 



3 



< lim sup ■ 



1 - <f)(e- x ) 

N 



F Ef=i(W(e- x A)) 



3 



1 - ^(e-*) 

< lim sup 3 



< lim sup ■ 



F 



(> + £f=i 



At t 



N 



Ov + Ef=i 



A, t 



iV + 5> 



1=1 



because 



+ Eili A)] < (Af + 1)3 [iV] < oo. 
Step 1, case ii) A; 's are unbounded. Take M > big enough that will be specified later. Define 
a truncated random vector (B(v),Ai(v), . . . A N (v)) = (B(v) AM, At(v) AM, . . . , A N (v) AM). Now 
take 



m M {t) = 3 



,i=i 



and observe that rriM(t) = 1 has two different solution <xm < a < Pm and both of them converge 
to a as M — > oo. Indeed, it follows from the fact that tom(ck) < 1 and for t G (0, a + S) different 
than a the Lebesgue Theorem gives that m,M(t) — > m(t) > 1. We will assume that /3m < a + 6. 
Define 



fl=53i( W ) J B(«), 



where L is defined in the same way as L but in terms of Ai(v). Clearly R < R a.s. Since M can 



be chosen in such a way that for any j the measure P log Aj £ du, Aj > 0, N > j 



has the same 
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support as P [log Aj £ du, Aj > 0, N > j] we can apply Lemma 1.7 a) for the random variable R and 
obtain F\R > t] ~ C+f~' 3ju for some positive C+. Hence, by Tauberian Theorem 



(3-12) 



l-0(f)~C + r(l-/3 M )^ M , 



where 0(f) = 0^(f) > 0(f). Therefore we can find C such that 1 - <j>(t) > Cot 13 *' for < f < 1. On 
the other hand, by Lemma 3.3 and Corollary 2.2, 1 — 0(f) < Cif 3 for f3 < a. This implies that for 
some C and any £ 



1 - 0(f) < C 1 - 0(f) 



/3/Au 



From (3.12) get 



>° l-0(f) 

hence for any < e < ad/2 we may find C\ such that 

1 - 0(st) 



1 - 0(f) 

for any f > and s > 1. Therefore, we have 



< C x s 



/ ~ \ /3//3m 

1 - 0(fs) < C ( 1 - 0(fs))J 

< C fCi (l - 0(f)) 

<c 5 



(l-£(f))"%<^>A 



-A/ 



Thus, since the function F is increasing, F(u)/u f* is bounded for /3m < 2/3 and 0(f) < 0(f) we 
have 



^(Ef =1 (i 



1 - (j}{e~ x ) 



Aj))) ^(Eti(l" H^Q- VAQ))) 



1 - (j){e- x ) 



< 



AT 



1 - 0(er 



i=i 



AT 



Now, applying Lemma 3.5 with = Af , r = and p = (one can choose sufficiently 

large M and /3m close to /3 such that the pair r, p satisfies asumptions of Lemma 3.5) we obtain 



3 



< C[ E[N 1+S ] +E 



N 



a+(5 
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Finally, by Lebesgue Theorem and since F(u)/u —> as u —> we deduce that 



3 



lim 



(e-«) 



= E 



F 



lim 



(Ef=iO 



(e-^i)) 



(e-») 



= 0. 



Step 2. Let us introduce a family of functions 



My) 



g(g + y) 



Dividing the equation 

by D(x), we obtain 
(3.13) 



D(x + y) = E [D(z + y + Y)} - G{x + y) 



K{y) = E [h x (y + Y)] - G Jf + V \ h x {y). 



D(x + y) 

Since D(y)e~ ay = 1 — <f)(e~ v ) is decreasing and D(y)e^ 1 ~ a ^ v = (1 — <fi(e~ v )) / e~ y is increasing the 
same holds for functions h x . Therefore we conclude that h x (y) < max{e Qy , e( a ~ 1 ' v } and that h x are 
equi-continuous on bounded sets. By Arzela-Ascoli theorem the set {h x } is relatively compact in 
the topology of uniform convergence on compact sets. Take now an accumulation point h as x — > oo 
i.e. h Xn — > h for some sequence x n — > oo. Passing to infinity with x n in (3.13) from Step 1 and the 
Lebesgue theorem we have 

h(y)=E[h(y + Y)}. 

Since any positive y-harmonic function is constant it yields that h(y) = h(0) = 1. Hence, h is the 
unique accumulation point and therefore D(x + y)/D(x) — > 1. 

Step 3. Finally we get rid of the additional assumption (3.7). We define B = min{l,i?} and 
consider 

R = L(v)B(v), 

i.e. R is defined in the same way as R in (1.2) but with B(v) replaced by B(v). Of course R < R, 
hence R is also finite a.s. and solves the equation 

N 

(3.14) X = d J2AiXi + B, 



with the vector (£?, j4i, A2, . . . ) satisfying hypotheses (1.10)-(1.14). 

Let 4> = cj)^ be the Laplace transform of R. Notice that (3.7) is satisfied for R and B replaced 
by R and B, i.e. 



1 



(3.15) 



-ti 



lim 

1 - <P(t) 



= 0. 



Indeed, observe first that E[R] =00, otherwise we would have 



E[R] = E 



N 



5> 



E[R] + E[B] 
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which is impossible since 1 < E 



Ai < oo. Thus in consequence 1 — cf>(t)/t tends to infinity 



as t goes to 0. Moreover, since B is bounded a.s., by the Lebesguc theorem and by Lemma 3.3 we 
have 



lim 



1 -E 




1 - 


4>(t) 



limE 



1 - e 



-IB 



t 



t 



E[B] ■ lim 



1 - ((>(t) <j)(t) 



0. 



which proves (3.15). Therefore, we may use the results proved in the first two steps of the proof 
saying that 

e «(a+y)(l _ 2( e -i x +y))) 

lim - = 1 



Hence 



e ax (l - 4>{e~ x )) 
1 - 4>{t) 



f" 



= L(l/t), 



for some slowly varying function L. 

Since R < R and <fr > 0, for < e < 5 we have 



lim -t— < lim = — = lim 



By Lemma 3.3 



Finally, we get 



s->o 1 - <p(t) f->o 1 _ *->o L(l/t) 

l-E[e- lB ] 



t a+e 



< C < 



1-E[e-* S 1 l-E[e- tB l t a+e 

lim ~~: — r — - = lim ■ 



t^O 1 - (j)(t) t-K) t a+e 1 - (j)(t) 

The last Proposition implies immediately the following results 
Corollary 3.16. Under assumptions of Theorem 1.9 we have the following: for any s > 



□ 



(3.17) 



hmI^ = S «. 



t->o 1 - <f)(t) 

In particular, the function L(t) = (1 — <j)(l/t))t a is slowly varying. 

Corollary 3.18. Under assumptions (1.10)-(1.14) we have E [R^] < oo for (3 < a and E [i?* 9 ] = oo 
for (3 > a. 

Proof. We have 

/•OO poo 

E [RP] =/3 t^- l ¥ [R> t]dt = C Q + (3 / t'^^L^dt 
Jo Ji 

what is finite if (3 < a and infinite if (3 > a. □ 
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3.2. Some properties of the function G. To prove our main results we will use the renewal 
theorem, therefore before we will proceed with the final arguments we have to prove some properties 
of the function G. 

Lemma 3.19. There exists e > such that the function e e ^G(x) is directly Riemman integrable. 
Proof. For any x and e < min{5, 1 — a} we have 



e ±EX \G(x)\ < e {a±e)x E 



N 



N 



^(l-cKe-^Vj-l + n^e-^ 



>±£)X E 



N 



l) IJ^e-M,-) 



eS a±£)x E 



N 



N 



^(l-0(e-^,))-l + n^( e ^ 



i=l 



e {a±e)x E 



N 



l-e- e 



i=l 



= fi(x) + f 2 (x), 
Notice first that since for 7 < a + e, by Lemma (3.3) we have 



e lx E 



N 



< e lx ( 1 — E 



-e~ x B 



\ < Cniin{e^,e^- a - i ) 1 } , 



the function /2 is directly Riemman integrable. 



Let us now examine the function f\. First, we will prove that f\ is integrable. For this purpose 
notice that if 7 > we have 



/ 



e^'E 



JV 



dx 



^(l-^e-^-l + n^e-^) 

i=l 

N 

^(l-0( e -^.))-l + e -£f =1 (i- 



< / e lx E 



C A,)) 



< E 



. Vi=i 



dx 



(e- x Ai)) dx 



The above expression is finite for 7 < a + e. Indeed, from the monotonicity of F on the positive 
half line and from Lemma 3.3, for /3 < a we can bound the last integral by 



3 



J e**F (°lt A ? E [ R0 ] dx 



= E 



N 



x / e"> x F{e- px )dx. 



To see that the expression above is finite we apply Lemma 3.5 with r = (a + S)/a and Xi = A r ^ . 
The second term is finite since F(t) < rnm{i, t 2 /2} and j3 < 7 < 2(3. Thus f\ is integrable. 
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Next we have to check that f\ satisfies (2.4). Inequality (3.10) implies that e ^ Q±e - )x /i(a;) is 
decreasing, hence also e~ x fi(x), since a + s < 1. Therefore, for any h > and x £ I n (h) we have 

fi(x)<e x - nh h(nh)<e h h(nh) 

and similarly we can estimate from below 

fi(x) > e x -( n+1 >h{{n + l)h) > e- h h((n + l)h). 

Since f\ is intcgrable, the scries 



J2h(nh) = Y,fi((n + l)h)<J2l I 



is finite. We can write 



</<> 



e h h{x)dx 



fi(x)dx < oo 



]T sup \f x {x)-f x ( y )\. h <Y,{z h h{nh)-e- h h{{n+l)h))-h 

neZ x,yei n (h) neZ 

< V (e fe - e" fe ) /i(n/i) • h < V (e h - e"' 1 ) / A^e^x 



= (e^-l) / h(x)dx. 

JR 

The last expression converges to as h goes to 0, thus f\ is directly Riemann integrable. 



□ 



Corollary 3.20. Functions xG(x) and G{x) are directly Riemman integrable. 

Proof. The corollary follows from the fact that (1 + |x|)|G(x)| < C (e~ £ ° x + e £ ° x ) \G(x)\, for some 
sufficiently large C. □ 

Corollary 3.21. If J G(x)dx = then the function G(x) — J_ G(y)dy is also dRi and satisfies 
jG(x)dx = -J xG(x)dx. 



Proof. For x < we have 

\G(x)\ < 



G(y)dy 



< / e £oX e- Eoy \G(y)\dy < Ce e 



In the same way we prove for x > that 



\G(x)\ < 



G(y)dy 



G(y)dy 



< Ce~ eoX . 



Hence G is directly Riemman integrable. Moreover, 



G(x)dx 



G{y)dydx 

r 

G{y)dydx 



x<0 J — oo 


G(y)dxdy 
V<0 Jy 



x>0 J x 



G(y)dydx 



y>0 JO 



G(y)dxdy = yG(y)dy. 



□ 
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3.3. Existence of the limit. We are able now to describe behaviour of <f> near 0. 
Proposition 3.22. Under hypotheses (1.10)-(1.14) we have 

limi^M = C + 

t^Q t a 

Thus, by the Tauberian theorem (Lemma 3.1) we deduce 
Corollary 3.23. Under hypotheses (1.10)-(1.14) we have 

lim x a P[R > x] = C+ 

x — >oo 

Proof of Proposition 3.22. The scheme of the proof is similar to the proof of Theorem 2.18 in [13]. 
Let {Y n } be a sequence of independent copies of Y defined in (2.5). By S n we denote their partial 
sums, i.e. S n = Y17=i We define the stopping time L = inf{?i > : S n < 0} and the sequence of 
stopping times = inf{n > T^-i ■ S n > S , T fc _ 1 }- 

Since D is a solution of the Poisson equation the sequence of random variables 

71-1 

M n (x) = D(x + S n )-J2 G(x + Si) 

forms a martingale with respect to the natural filtration generated by {5^}- In view of the optional 
stopping theorem we have 

E[M nAL {x)]=E[M {x)]=D(x) 

and equivalently 

~nAL— l 



E[D(x + S nAL )]-D(x) =E 



G(x + Si) 



i=0 



Next we want to pass with n to infinity. Notice that by the duality principle [14] 



E 



L-l 



Y,\G(x + Si)\ 



E 



Y,\G(x + S Ti )\ 



and the last sum is finite since G is dRi, hence we can pass with n — > oo on the right side. In order 
to justify passing to limit on the left side observe that by Proposition 3.6 for any e < S we have 
D(x) < Ce^L Since E [e eSi ] < oo (see [14] (3.6a) in Chap. XII), we can pass to infinity. 

Thus we obtain 

L-l 

Y,G(x + S l ) =:R(x). 



i=0 



(3.24) E[D(x + S L )]-D(x) = E 
Applying again the duality principle we have 

oo 

(3.25) R(x) = ^G(x + S T J- 
Now the renewal theorem yields that 

(3.26) lim R(x) 

x— foo 

Integrating (3.24) we have 



n=0 



J M G(x)dx 



(3.27) 



[S Tl ] 

E[D(y + S L )]-D(y))dy = 



R{y)dy 
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what can be rewritten as 
(3.28) D(x) ■ E 



D{x + y) 
D(x) 



dy 









- E 


/ " D{y)dy 


-I 




Jo 





R(y)dy 



By Proposition 3.6, D{x + y)/D(x) < Ce ey hence again the same argument as before tells us that 
we can pass to the limit under the integral sign and obtain 



lim E 

x— >oo 



D(x + y) 
D(x) 



dy 



E[S L 



Dividing by x in (3.28) and passing to the limit we obtain 



and finally 
(3.29) 



D(x) r , LG(x)dx 
lim — E [Sl] = lim R(x) = - \,y, , 



lim 



a; — y oo x 



E[S Ti ]E[Sl] 



where cr 2 := Vary = 2E [-S L ] E [S Tl ] (see the proof of T18.1 on page 196 in [25] 
Notice that in view of Lemma 3.1 the last formula implies 



lim 



t" 



»[R > t] 



2j R G(x)dx 



t-)-cx) logt " 

By Proposition 2.1 the last constant must be 0, therefore 



G(x)dx = 0. 

Now we repeat the above procedure. Integrating (3.24) we have 



E[D(y + S L )]-D(y))dy 



E 



i=0 



dy 



what is equivalent to 
E 



D(x + y)dy 
Passing with x to infinity in 

D(x)E 
we obtain 

lim D(x) 



oo „ x 

G(y + S Ti )dy 



J2 G{x + S Tt 



i=0 



D(x + y) 



dy 



J2g{x + s Tz ) 



i=0 



D(x) 

_ 2 jG(x)dx -2 J xG{x)dx 



□ 
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3.4. Positivity of the limiting constant. Now we proceed with the last step of the proof and we 
justify that the constant C+ in the statement of Corollary 3.23 is strictly positive. We follow here 
the method developed in [11]. 

The key tool to deal with this problem is the Landau Theorem. Originally it was stated for 
Dirichlet series but it can be extended for functions of the type f(s) = J x s /i(dx) where \i is a 
positive measure on R + (compare with [27, Theorems 5a and 5b in Chap. II]) 

Theorem 3.30 (Landau). Let a c be the abscissa of convergence for f{z) = J °° x z [i{dx\ i.e. the 
integral converges for 5Rz < er c and diverges for 5Rz > a c . Then o c is a singularity for f . 

As a conclusion we get that if z M> E [R z ] has an analytic extension on some open neighborhood 
of a then E [R z ] is well defined there. Our aim is to find such an extension, under the assumptions 
that / G(x)dx = and / xG{x)dx = 0. 

Theorem 3.31. Under assumptions (1.10)-(1.14) the constant C+ is strictly positive. 
Proof. Let us define 



By Lemma 3.19, H is well defined for < 3?z < a + eq. Moreover, H is a holomorphic function 
for < < a + e . Indeed, take any closed piecewise C 1 curve 7 in < < a + e . Then by 
Lemma 3.19 we have 



H(z)= / e xz E ^(1 - <j>(e- x Ai)) - (1 - ^(e" x )) da. 






Hence, by the Fubini Theorem 




Therefore, by Morera Theorem H is holomorphic in < 3ft z < a + Eq and 




where the second equality follows from (e 



,xh 



l)/h < Ce e ' x ' and the Lebesgue Theorem. 
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Now take z such that Kz < a and observe that for any a > 

"aR 

e zx (l - (j){e~ x a))dx 



E 



L'" K \L 



l du 



dx 



Je x <aR/u 

a 2 E [R z }T(l-z" 



= E 


~a z R z r 






. z Jo 



e u u z du 



< 00. 



Therefore for Kz < a we can express H in another form: 



H(z) 



' N 

.i=i 



(W( 



(Z.r 



e"*^)) - (1 - ^(e-)) 

- 1) r(i - z) E [jj.] ( TO ( 2 ) _ X ) r( X _ z) 



From this we obtain that for < iffiz < a 
(3.32) E [i? z ] 



zH(z) 



m{z) - l)r(l - z) 

However, if we assume that both H(a) = J G(x)dx and H'(a) = / xG{x)dx are equal to 0, then a 
is a root of multiplicity at least two of H. By our assumptions m(a) — 1=0, m'(a) = 0. The strong 
convexity of m yields m"(a) > 0. Since T(l — a) ^ the right side of (3.32) defines a holomorphic 
function on < < a + eo that extends E [R z ]. From the Landau Theorem we know that if E [R z ] 
has a holomorphic extension on 5Rz < a + Eq then it is expressed by the same formula E [R z ]. In 



particular, one can find positive 5' such that E 



< oo, contrary to the Corollary 3.18. 



□ 
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